Mon. Not. R. Astron. Soc. 000, [|-| (2001) Printed 1 February 2008 (MN WI^ style file vl.4) 



The effects of seeing on Sersic profiles. II. The Moffat PSF 



I. Tmjillo,*^ J. A. L. Aguerri,^ J. Cepa^ and C. M. Gutierrez 

^ Instituto de Astrofisica de Canarias, E-38205 La Laguna, Tenerife, Spain 

^ Astronomisches Institut der Universitdt Basel, Venusstrasse 7, CH-4102 Binningen, Switzerland 



Accepted 0000 December 00. Received 0000 December 00; in original form 0000 October 00 



o 
o 

(N 

CD 

in 



> 

\o 
o 
o\ 
o 



ABSTRACT 

The effects of seeing on Sersic r^/" profile parameters are extensively studied using 
a Moffat function. This analytical approximation to the point spread function (PSF) is 
shown to provide the best fit to the PSF predicted from atmospheric turbulence theory 
when (3 ^ 4.765. The Moffat PSF is additionally shown to contain the Gaussian PSF 
as a limiting case (/? — > oo). The Moffat function is also shown to be numerically well 
behaved when modelling narrow PSFs in iKT images. Seeing effects are computed for 
elliptically symmetric surface brightness distributions. The widely used assumption of 
circular symmetry when studying the effects of seeing on intrinsically elliptical sources 
is shown to produce significant discrepancies with respect to the true effects of seeing 
on these sources. A prescription to correct raw (observed) central intensities, effective 
radii, index n and mean effective surface brightness is given. 

Key words: atmospheric effects - methods: data analysis - galaxies: distances and 
redshifts - galaxies: photometry. 
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1 INTRODUCTION 

It is well known that the ability to parameterize galax- 
ies from ground-based images is severely compromised by 
seeing, which scatters light from the objects producing a 
loss of resolution in the images, lower mean surface bright- 
nesses than the true values and larger effective radii. The 
effects of seeing have been extensively studied in the case 
of elliptical galaxies with r^^* profiles (Franx, Ilhngworth & 
Heckman 1989; Saglia et al. 1993). Recently, Trujillo et al. 
(2001, hereafter TOl) extended previous work by studying 
analytically the effects of seeing on elliptically symmetric 
surface brightness distributions following the Sersic (1968) 
r^^" law and assuming a Gaussian point spread function 
(PSF). Sersic's generalization of the de Vaucouleurs (1948, 
1959) r^/* law has been shown to provide a better represen- 
tation to the distribution of light in both elliptical galax- 
ies (including the dwarf ellipticals) and the bulges of Spiral 
galaxies (Caon, Cappacioli, & D'Onofrio 1993; D'Onofrio, 
Capaccioli, & Caon 1994; Young & Currie 1994; Andredakis, 
Peletier, & Balcells 1995). 

The existence of "wings" in stellar profiles reveals that 
the real PSF deviates from the Gaussian form. In this paper 
we show, from the size of the wings present in real images 
(e.g. Saglia et al. 1993), that such deviations from Gaussian 
PSFs can result in difi'erent values for the profile parame- 
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ters in the range of 10-30%. The new generation of ground- 
based telescopes and the study of galaxies at high redshifts 
make these types of studies crucial in order to obtain reli- 
able (unbiased) information from the structural analysis of 
these objects. 

This paper presents a further, more detailed, analysis 
of the effects of seeing on Sersic profiles when "wings" are 
present in the PSF. For this reason we have modelled the 
PSF by a generalization of the Gaussian form: the Moffat 
function (Moffat 1969), which describes well the presence 
of wings. It should be noted that these kinds of studies are 
not only important for ground-based observations. In fact, 
HST images present their own "narrow" PSFs (see a de- 
tailed study in Bendinelli, Zavatti & Parmeggiani 1987, and 
Krist 1993). The use of these steep PSFs presents numeri- 
cal problems which can be avoided by modelling the narrow 
PSFs with polynomials instead of exponential expressions 
like Gaussians. In Section 2 we summarize some general re- 
sults from the use of Moffat PSFs. Section 3 describes the 
effects of seeing on the Sersic profile parameters brought 
about by the Moffat PSF. A prescription for seeing correc- 
tions is given in Section 4. 



2 GENERAL REMARKS ABOUT MOFFAT 
CONVOLUTION 

Point spread functions can be determined observationally 
by studying the scattering of stellar light. Numerous papers 
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have been devoted to this problem (e.g. Moffat 1969; King 
1971; Bendinelh et al. 1990). Among the analytical approx- 
imations, the Moffat function (see Eq. [1]) has been widely 
used to model the PSF (e.g. Bendinelli, Zavatti & Parmeg- 
giani 1988a, b; Young et al. 1998); for instance, the IRAF 
data reduction package (Tody 1986) adopts the Moffat func- 
tion as a standard PSF. In Fig. 1 we plot both the Gaussian 
function and Moffat function having a range of /3. Note that 
as 13 increases the Moffat function tends to approximate the 
core of the Gaussian profile. In fact, a Moffat function con- 
tains the Gaussian PSF as a limiting case (see Appendix A) . 
Moreover, the Moffat PSF has two clear advantages over the 
Gaussian PSF: 



elliptical source depends only on ^: 7(r) = /(^). The con- 
volution equation that represents the effect of seeing on the 
surface brightness distribution is given by: 

Jo Jo 
where PSF(c',e»',,t,6») is the Moffat PSF given by: 

PSF(c',e',e,^) = ^ 1 + 

+ - 2^C'cos(6l - e') + (e^ - 2e)(C' smO' - ^ sin 61)2" 



(i) It is numerically well behaved in the treatment of nar- 
row PSFs. 

(ii) It allows the "wings" that usually appear in stellar 
profiles to be fitted. 

Very accurate convolutions between the PSF and the 
model profiles of the galaxies are required in order to ob- 
tain reliable results. Current reduction packages use Fast 
Fourier Transforms to evaluate the convolutions. This is in 
fact inappropriate where there are strong changes in the in- 
tensity gradients of the galaxy profiles. The inner parts of 
galaxy profiles are steep, and this demands a very accurate 
measurement of the high frequencies in the Fourier domain. 
Narrow PSFs (such as those of the HST) magnify this prob- 
lem because they also present a steeper profile. Working in 
the real domain does not exempt us from trouble either; in 
fact, one can encounter several numerical problems when 
performing accurate convolutions when using a Gaussian to 
model narrow PSFs. Current computers can manage num- 
bers of the order of ~ e^"". These kinds of numbers can be 
eeisily obtained when working with Gaussians which have 
cr < 1 in units of pixels (see the exponential expressions at 
play when performing a Gaussian convolution in Eq. [4] in 
TOl). The use of Moffat functions avoids this problem due 
to the use of polynomials instead of exponential expressions 
(see Eq. [3]). In this sense, Moffat functions are numerically 
better behaved than Gaussians when dealing with narrow 
PSFs. 



2.1 Mathematical analysis 

We will use a circular Moffat function to model the point 
spread function: 



PSF(r) = 



1 + 



(9' 



(1) 



with the fuU width at half maximum, FWHM=2aV2i/'' - 1, 
where PSF(FWHM/2) = (1/2)PSF(0) and the total flux is 
normalized to 1. Consider a case where, in the absence of 
seeing, the surfaces brightness distribution, /(r), of a galaxy 
is elliptically symmetric. This means that the isophotes of 
the object all have the same constant ellipticity e (e = 1 — 
b/a, where a and b are respectively the semi-major and semi- 
minor axes of the isophotes). 

As shown in TOl elliptical coordinates (^, 6) are the 
most appropriate for this type of problem. In this coordi- 
nate system, the surface brightness distribution, /(r), of an 



The subscript "c" shall be used from here on to re- 
fer seeing-convolved quantities. Along the major axis of the 
object, 6=Q, the angular integral can be solved analytically 
(6>0): 



dO PSF(C ,61 ,C,0) = 2 



/3-1 
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k = 

n 2k + l\ 
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where 



= 1 H and w 



1 C 1 

(2e- 62)1/2 0^1 



/2 



(4) 



(5) 



and Cn{t) and B{z,w) are the Gegenbauer polynomials 
(Gradshteyn & Ryzhik 1980, p. 1029) and beta functions 
(Abramowitz & Stegun 1964, p. 258) respectively. A simpler 
expression is obtained for the convolution in the circularly 
symmetric case (e = 0): 



7c(r) = 2^-^/ dr'r'l{r')ao^Pfi.i 
Jo 

where 



1 / , ^ r^ + r^ 
an IH ^ — 



ao 



/2 „ ^2 



+ 



1/2 



(7) 



and Pn{x) is the Legendre function of first class 
(Abramowitz & Stegun 1964, p. 332). 

2.2 The effect of seeing on the central intensity 

For any intensity distribution with elliptical symmetry, /(^), 
the seeing convolved central intensity, Ic(C = 0), is such that 



/c(0) 



where 



fe=l + ^ 



(1-e) / e /(c ) 



and c ; 



(62 _ c2)'9/2 



/3-1 



-^(^^-2.). (9) 



Note that the effect on the central intensity of the seeing 
convolved distribution is a function of the intrinsic ellipticity 
of the object. Basically, as e increases the spread of photons 
from the inner parts of the profile due to the seeing is more 
efficient and consequently the central intensity decreases. 
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Figure 1. Top panel: The normalized intensity profile of a Gaussian and Moff^at PSF function having different values of (3 are plotted 
against the radius of the PSF in units of FWHM. Bottom panel: The difference between the normalized Moffat PSF (PSF^(r)) and the 
normalized Gaussian PSF (PSF(j(r)). 



2.3 The effect of seeing on the ellipticity of the 
isophotes 

In the absence of seeing, by construction, all isophotes of 
the profile have the same ellipticity, whereas the presence 
of seeing tends to make them circular. Using the isophote 
condition, /c(^,0) — Ic{^,n/2) it is possible to derive an 
implicit equation that gives the variation of the ellipticity 
with the radial distance: 



1 2k + 1 



B 



0. 



(10) 



For this problem, it is useful to introduce a^,w,a* and 
w* as functions of the Cartesian coordinates x, y: 



= 1 H and w 



and 

= i H and w 



1 X 1 

\l/2 r. 1/2 



(2e - £2)1/2 a ^1/2 



'(2e- £2)1/2 Q, (a|)i/2- 



(12) 



From this implicit equation we can obtain y/x and therefore 
the ellipticity of the isophotes affected by seeing using e{x) = 
1 - y/x. 



kJo(kr) 



<^-rwim-) dk. 



(13) 



where Jo is the standard Bessel function (Abramowitz & 
Stegun 1964, p. 358). For a given FWHM, we have evalu- 
ated the value of /3 that minimizes the difference between 
the prediction of the atmospheric turbulence theory and the 
Moffat function by minimizing the of the fit between both 
PSFs. An optimum value of /3 ~ 4.765 was found. In Fig. 
2 we have shown the difference between the PSF prediction 
from turbulence theory and a Moffat function for a value of 
P — 4.765. It can be seen that the agreement is quite good. A 
Moffat function could therefore be used to reliably model the 
turbulence prediction, although, the PSFs usually measured 
in real images have bigger "wings", or equivalently smaller 
values of /3, than those expected from the turbulence theory 
(e.g. Saglia et al. 1993). This is because the real seeing not 
only depends on atmospheric conditions but is also caused 
by imperfections in telescope optics. 

The presence of these bigger "wings" in real images 
makes Moffat functions a better choice to model the PSF 
than the turbulence theory prediction. As an example of 
this, current packages of data reduction in IRAF suggest a 
default value of /3 = 2.5. In order to span the range of the 
different "wing" sizes present in real images, in the next sec- 
tion we have modelled the Moffat PSFs using three different 
values of /3: p—5 (to simulate the turbulence prediction), 2.5 
(the default value of the IRAF package), and 1.5 (to model 
a large "wing" in the PSF). 



2.4 The abihty of the Moffat function to match 
the atmospheric turbulence prediction of the 
PSF 

The theory of atmospheric turbulence predicts the PSF to 
be the Fourier transform of exp [-(kbf^^] (Fried 1966; Woolf 
1982), where FWHM = 2.9207006& and 6 is a scaling param- 
eter. In the real domain this PSF is written as: 



3 THE EFFECTS OF SEEING ON THE SERSIC 
PROFILE PARAMETERS 

The equations that we have shown in the previous section 
are general results for Moffat seeing. For practical purposes 
with applications to real galaxies, we are going to focus on 
the Sersic profile. In the particular case of r^^" profiles, the 
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Figure 2. Top panel: The best fit to the PSF predicted by the theory of atmospheric turbulence (diamonds) using a Moffat function 
with a /3 value of 4.765 (solid line). Bottom panel: The difference between these two PSFs as a function of the radius of the PSF in unit 
of FWHM. 



surface brightness distribution is given (in elliptical coordi- 
nates) by: 

/(C) = /(0)10-'"[<^'''''=)*i, (14) 

where /(O) is the central intensity and the effective radius 
of the profile. The constant &„ is chosen such that half the 
total luminosity predicted by the law comes from ^ < re- bn 
can be well approximated by the relation b„ — 0.868n— 0.142 
(Caon et al. 1993). 

3.1 The central intensity 

To study the effect of seeing on the central intensity we 
make use of Eqs. (8) and (14). Figure 3 shows this effect for 
different values of the ellipticity in the intrinsic light profile, 
assuming /3 = 1.5, 2.5, 5 and for the Gaussian case (i.e. 
/3 — > oo) as well. 

na As P increases (i.e. as the size of the PSF "wings" 
decrease) we asymptotically recover the effect on the central 
intensity produced by a Gaussian PSF. 

From this figure it follows that the effect of seeing on 
the central intensity increases as the size of the "wings" in- 
crease. This is easily understood as broader "wings" increase 
the probability that a photon will hit the imaging device at 
a point further offset from where it would have hit in the 
absence of a seeing. For a typical value of /3 (e.g. /3 — 2.5), 
the difference from a Gaussian PSF is ~ 10%. 

For a given seeing-disc (n.b. the relative size of the 
FWHM to effective radius increases along the x-axis of Fig. 
3, c.f. Fig. 4 and Fig. 5), the central intensity of profiles with 
larger values of n is more affected than for low n as is ex- 
pected because of the higher central hght concentration of 
these profiles. As noted in Section 2.2, the effect of seeing 
on the central intensity of the object is also dependent on 
the intrinsic ellipticity of the object: the central intensity of 
galaxies with larger ellipticities are more affected by seeing. 



3.2 The effective radius 

The seeing effect on effective radius can be obtained by solv- 
ing for rj; from the conservation of luminosity by the con- 
volution L'^{r^)=L{rs), where L{re) is the luminosity of the 
source inside and L'^ir^) is the luminosity obtained from 
the object affected by seeing, measured inside its effective 
radius. 

Figure 4 shows this effect for different values of elliptic- 
ity, with P = 1.5, 2.5, 5 and for the Gaussian PSF. Here, 
the presence of significant "wings" in the PSF produces an 
effective radius larger than what is expected from Gaussian 
seeing. As n increases the convolved effective radius also in- 
creases. The ellipticity effect is also shown. Greater elliptic- 
ities result in greater effective radii, and these differences 
are more important for greater values of n. This result is 
as expected due to the diminution of the central intensity 
with larger ellipticity. For the values of P typically present 
in real images (2.5 < /3 < 4; see Saglia et al. 1993) we obtain 
deviations from Gaussian seeing in the range 15-30% (big- 
ger deviations are obtained for smaller values of the ratio 
r=/FWHM). 

ft should be noted that our measurement of the effective 
radius has been obtained over the semi-major axis. Some 
authors use as radial distance the magnitude r* — ^fab; in 
this case, the effective radius of the object affected by seeing 
is given by r" — r^^Jl ~ ^ ('"c), where e{r'^) can be obtained 
using Equation (10). 



3.3 The Sersic index n 

To quantify the effect of seeing on the shape parameter n we 
use the parameter ?7(C) (TOl). r){X) is equivalent, locally, to 
the parameter n of the Sersic profile. This parameter can be 
understood as a measure of the slope of the profile. In Fig. 
5 we show the effects of seeing on this parameter (evaluated 
at rj;) for different sizes of the PSF "wings" and different 
intrinsic ellipticities. It is easy to see how the real value of 
n is recovered asymptotically when the ratio rJ;/FWHM in- 
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Figure 3. The effects of seeing on the observed central intensity, /c(0), for different values of n. The Gaussian case (solid line) and three 
different values of /3 are shown: (3 = 5 (dot-dashed line), (3 = 2.5 (dashed line) and /? = 1.5 (dot-dot-dot-dashed line). Three different 
ellipticities are also shown: e = (left column), e = 0.25 (middle column) and e = 0.5 (right column). 



Figure 4. The effects of seeing on the observed effective radius, for different values of n. The Gaussian case (solid line) and three 
different values of /? are shown: /3 = 5 (dashed line), /3 = 2.5 (dotted line) and /3 = 1.5 (dot-dashed line). Three different ellipticities for 
the sources are also shown, e = (left column), e = 0.25 (middle column) and e = 0.5 (right column). 



creases. The effect of seeing on ri{r'^) is bigger for smaller 
values of l3. This means that bigger PSF "wings" produce 
a stronger effect on the slope of the profile. The increase 
in the intrinsic ellipticity of the profile has a similar effect, 
but the influence is not as important as it was for the pre- 
vious parameters. Note that seeing effects always produce 
a surface brightness profile with a smaller value of n than 
the actual one. It is expected that any procedure to recover 
the structural parameters of the profile that does not take 
into account the effect of seeing will obtain lower values of n 
than the actual ones. Lower values of n will also be expected 
if the intrinsic ellipticities of the objects are not taken into 
account during the recovery process. This is crucial in the 
study of high-z galaxies. As was shown for the central in- 
tensity and effective radius, the presence of "wings" in the 
PSF causes the values of the profile parameters to deviate 



from the prediction made correcting for Gaussian seeing. In 
the case of the ri{rl) parameter, these deviations are in the 
range of 10-20%. 

3.4 The mean effective surface brightness 

From the previous results it is now easy to study the effect 
of seeing on the mean effective surface brightness, defined 
as; 

(M)e = -2.51og%^. (15) 

Try c 

By using Eq. (15) and the conservation of the flux, it imme- 
diately follows that 

A{^*)e = (M)c-{M)c = 51og^. (16) 
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Figure 5. The effects of seeing on the index for different values of n. The Gaussian case (solid line) and three different values of /3 

are shown: /3 = 5 (dashed line), fi = 2.5 (dotted line) and (3 = 1.5 (dot-dashed line). Three different ellipticities for the sources are also 
shown, e = (left column), e = 0.25 (middle column) and t = 0.5 (right column). 



Figure 6. The differences A(^)o between the measured mean effective surface brightness, and the seeing-free quantity (/i)o for 

different values of n. The Gaussian case (solid line) and three different values of /3 are shown: (3 = 5 (dashed line), /? = 2.5 (dotted line) 
and 13 = 1.5 (dot-dashed line). Three different ellipticities for the sources are also shown, e = (left column), e = 0.25 (middle column) 
and e = 0.5 (right column). 



In Fig. 6 we show how the mean eflFective surface bright- 
ness changes as a function of r^/FWHM for different val- 
ues of n and intrinsic ellipticities. This figure clearly shows 
that galaxies affected by seeing have apparent mean surfaces 
brightnesses lower than their true values. Lower values of /3 
produce greater effects on this quantity. Also, as the intrin- 
sic ellipticity of the object increases the effects of the seeing 
on the mean effective surface brightness also increase. 



4 A PRESCRIPTION FOR SEEING 
CORRECTIONS 

It is possible to obtain the parameters of the Sersic profiles 
(seeing-free quantities) from the convolved quantities. We 



present an easy prescriptionf] based on the use of the plots 
of Figures 3, 4 and 5. The steps which observer must follows 
are: 

• Determine the FWHM and the value of /3 from the stel- 
lar profile by fitting a Moffat function. Current astronomical 
data reduction packages, such as IRAF, allow this fitting. 

• Measure rj; along the semi-major axis directly from the 
raw images. This can be done by solving the implicit equa- 
tion L^irl) = (l/2)L'^(oo). 

• Determine ?7(rJ) numerically using the expression: 



t A similar prescription for Gaussian seeing is presented in TOl. 
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1 I.iK) Jc(rg) 
J.(0) 



(17) 



• Evaluate the value of n and e from the use of Figure 5. 

• Obtain the value of using Figure 4. 

• Obtain the value of 7(0) using Figure 3. 



5 CONCLUSIONS 

As redshift increases, the apparent size of galaxies becomes 
progressively smaller and the affect of seeing progressively 
stronger. As we have seen, precise corrections for seeing are 
demanded in order to obtain reliable and comparable in- 
formation about structural parameters from objects at the 
same or different redshifts. We have chosen the Moffat func- 
tion to model the PSFs of real images. This choice has been 
made by following two criteria: the ability of this function 
to model the "wings" of the PSFs present in real images ob- 
tained from ground-based telescopes, and it is well behaved 
numerically because of its polynomial structure. 

We have studied the general properties of the Moffat 
function when modelling PSFs. These properties can be 
summarized as follows: it is a very good option to model 
the narrow PSFs present, for example, in HST images be- 
cause it is numerically well behaved; the Gaussian PSF is a 
limiting case of the Moffat PSF (/? — + oo); and the prediction 
for the PSF due to the theory of atmospheric turbulence can 
be numerically well approximated by a Moffat function with 
/3 ~ 4.765. 

For practical purposes we have analysed the effects of 
seeing caused by this PSF on the Sersic model. The effects 
on the central intensity, effective radius, n index and mean 
effective surface brightness are extensively shown in Figures 
3, 4, 5 and 6. We have also given an easy prescription for 
seeing correction that can be useful for observers in order to 
obtain the seeing- free quantities. 

Our main results have been to show the importance of 
taking into account the intrinsic ellipticities of the objects 
and the presence of "wings" in the PSFs for the recovery 
of accurate structural parameter. It is not sufficient to con- 
sider the PSF as Gaussian and assume circular symmetry to 
model the effects of seeing on the surface brightness distri- 
bution when the ratio of the effective radius to the FWHM 
is small (< 2.5). 
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APPENDIX A: THE GAUSSIAN AS A 
LIMITING CASE OF THE MOFFAT FUNCTION 

The Gaussian function can be obtained from the Moffat 
function (MF) in the limiting case where (3 ^ oo. One can 
rewrite the MF as a function of FWHM (F) and (3: 

-13 



PSF(r) =4(2 



1//3 , 



rf2 



1+4(2 



1//3 



.(Al) 



As /3 — > oo, we can substitute 2^^^^ — 1 with (ln2)//3, so: 

oi -0 



hm PSF(r) = hm — — - 

/3-.CX) 0^oo (3 TtF^ 



1 + 



4 In 2 



P 



Using limm-.oo(l + z/m)"^ = we have 

hm PSF(r) = —e . 

Finally, writing = 8a^ In 2 we obtain: 



lim PSF(r) = —— 



e 2 



.(A2) 



(A3) 



(A4) 



For practical purposes, a value of f3 — 100 is completely sat- 
isfactory for modelling a Gaussian by using a Moffat func- 
tion. 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science M^SX style file. 
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